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Abstract 

Periodic oscillations play an important role in many biomedical systems. Proper functioning of biological systems that 
respond to periodic signals requires the ability to synchronize with the periodic excitation. For example, the sleep/wake 
cycle is a manifestation of an internal timing system that synchronizes to the solar day. In the terminology of systems 
theory, the biological system must entrain or phase-lock to the periodic excitation. Entrainment is also important in 
synthetic biology. For example, connecting several artificial biological systems that entrain to a common clock may lead to a 
well-functioning modular system. The cell-cycle is a periodic program that regulates DNA synthesis and cell division. Recent 
biological studies suggest that cell-cycle related genes entrain to this periodic program at the gene translation level, leading 
to periodically-varying protein levels of these genes. The ribosome flow model (RFM) is a deterministic model obtained via a 
mean-field approximation of a stochastic model from statistical physics that has been used to model numerous processes 
including ribosome flow along the mRNA. Here we analyze the RFM under the assumption that the initiation and/or 
transition rates vary periodically with a common period T. We show that the ribosome distribution profile in the RFM 
entrains to this periodic excitation. In particular, the protein synthesis pattern converges to a unique periodic solution with 
period T. To the best of our knowledge, this is the first proof of entrainment in a mathematical model for translation that 
encapsulates aspects such as initiation and termination rates, ribosomal movement and interactions, and non- 
homogeneous elongation speeds along the mRNA. Our results support the conjecture that periodic oscillations in tRNA 
levels and other factors related to the translation process can induce periodic oscillations in protein levels, and may suggest 
a new approach for re-engineering genetic systems to obtain a desired, periodic, protein synthesis rate. 
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Introduction 

External and internal periodic oscillations play an important 
role in intracellular and extracellular biomedical systems and have 
attracted enormous research interest (see e.g. [1] and the 
references therein). Proper functioning of cells that are exposed 
to such periodic signals requires internal biological mechanisms 
that are able to synchronize with the periodic excitation. In the 
terminology of systems theory, the biological system must entrain or 
phase-lock to the periodic excitation. In other words, in response 
to a periodic excitation with period T the system's internal state 
converges to a periodic signal with period T. 

Entrainment in biological systems (sometimes called phase 
locking [2]) and, more generally, biological oscillators and rhythms 
have recently attracted enormous attention (see e.g. [3-5] and the 
references therein). For example, the sleep/wake cycle is a 
manifestation of an internal timing system that entrains to the 
24 hours period of the solar day using a visual pathway connecting 
the retina to the suprachiasmatic nucleus (SCN) [6]. 

Entrainment is also important in synthetic biology. For example, 
most hormones in the body are released in periodic pulses. 
Glucocorticoid secretion, for instance, has a circadian and 



ultradian pattern of release. Synthetic biological oscillators may 
be used to mimic these periodic release patterns in the 
administration of synthetic hormones to patients suffering from 
glucocorticoid-responsive diseases, thus improving therapeutic 
effectiveness [7]. 

The design of robust synthetic biological oscillators is also the 
first step for applications such as clocks that synchronize in order 
to coordinate intracellular behavior, and artificial platforms that 
can measure the genomic response to an oscillatory excitation [8] . 

Entrainment at the intra-cellular gene expression level 

Proteins are "tiny machines" performing a vast array of 
functions within living organisms including intra- and inter- 
cellular signaling transduction, immunological response against 
pathogens, movement of cells and tissues, facilitation of biochem- 
ical reactions, structure and support of the cell and tissue, and 
transport. Regions in the DNA, called genes, encode the 
information needed to produce proteins. Gene expression is the 
process by which the information inscribed in the genes is 
converted into proteins. The major steps of gene expression are 
transcription, translation, and mRNA and protein turnover [9]. 
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The cell-cycle is a periodic program that regulates DNA 
synthesis and cell division. Proper execution of the cell-cycle 
requires the expression and activation of key proteins at specific 

times along the period. This process must be tightly regulated, as 
perturbations in cell-cycle progression can lead to apoptosis or 
cancer. 

Recently, there is growing evidence that protein levels of cell- 
cycle related genes can be regulated not only via transcription (see, 
for example, [10]) but also via the translation elongation step. 
Higareda-Mendoza and Pardo-Galvan [1 1] investigated the role 
of human translation initiation factor 3 (eIF3) in cell-cycle control 
of A549 cells. They reported that eIF3f expression oscillates during 
cell-cycle, with one maximum expression peak in the early S phase 
and a second during mitosis. Their interpretation is that eIF3f is a 
translational modulator that selects mRNAs at specific cell-cycle 
phase time points. 

Frenkel-Morgenstern et al. [12] have shown that cell-cycle 
regulated genes tend to include non-optimal codons, i.e., codons 
that arc rare and are recognized by tRNA molecules with low 
intra-cellular abundance, and thus with a low translation rate. 
These codons create "botdenecks" in the translation process and 
thus their slow translation rate becomes rate limiting. They argue 
that periodicity in the tRNA levels of these codons induces 
periodicity in the translation rate of these genes. The fact that cell- 
cycle regulated genes display different codon preferences at 
different phases of the cell-cycle supports the conjecture that cells 
exploit non-optimal codons to generate cell-cycle-dependent 
dynamics of proteins via the translation process. In other words, 
the translation process entrains to the excitation generated by 
periodically var^'ing tRNA levels. 

In another recent study, Patil et al. [13] ha\'c reported an 
additional mechanisms by which cell-cycle can be regulated via 
translation. The ribonucleotide reductase (RNR) complex plays an 
important role in regulating cell-cycle transitions and in DNA 
damage response. They showed that the levels of 16 tRNA 
modifications, and thus the translation efficiencies of different 
codons oscillate during cell-cycle; in addition, these oscillations match 
the wobble interaction needed for translating codon of genes such 
as RNRl. Their results imply that translation regulation has a 
direct role in cell-cycle related oscillations. 

Two other recent studies [14,15] suggest that non-optimal 
codon usage during translation affect the expression, structure, 
and functioning of proteins, and are particulary important in the 
context of circadian clocks. 

Periodicity in gene expression that is related to periodic 
processes, such as the cell-cycle and biological clocks, is regulated 
at all the different gene expression stages. This includes 
transcription, translation, and post-translational regulation. The 
related regulation mechanisms include dozens (or even hundreds) 
of genes and proteins that interact with each other in ways that we 
are only beginning to unveil. Usually these networks of interac- 
tions include a few redundant mechanisms of oscillation regulation 
[14-20]. For example, it was suggested that cell-cycle regulation 
includes negative feedback oscillators. These can include for 
example the interconnection of two genes where the first gene up 
regulates the second, and the second down regulates the first [16]. 
Another possible regulation mechanism is via control of the 
transcription rate of tRNA genes (and other genes), resulting in 
oscillations in intra-cellular tRNA levels [12]. Since the decoding 
time of codons is affected by the available levels of the tRNA 
molecules recognizing them (see, for example, [21,22]), this may 
eventually lead to oscillations in the decoding times of different 
codons. 



It may sc-em natural to assume that periodic variations, with 
period T, in the initiation rate and/ or the decoding times of 
difierent codons will lead to a periodic protein production rate 
with the same period T. However, this assumption is actually quite 
strong. Indeed, these factors affect the protein synthesis rate via the 
dynamics of the translation mechanism, and not every dynamical 
system entrains to periodic excitations. Here we analyze this 
problem using a computational model for translation-elongation. 

Entrainment in a computational model for translation 

High-throughput experiments provide more and more data on 
the translation process. Computational models of translation are 
needed to organize, understand, and connect this data to various 

biophysical aspects of translation [23-28]. Understanding the 
dynamics of gc'iu" cxpri^ssion, and not only the static information 
encoded in the genes, is vital in order to und(;rstand how the 
biological components work together to comprise functioning cells 
and organisms. Developing a deeper understating of the dynamics 
of translation may thus have implications in many fields of science 
including human health [29-35], biotechnology [36^0], evolu- 
tion [31,41-47], functional genomics [48-54], and sy.stems biology 
[33,51,55-59]. Recent reviews related to translation may be found 
in [38,42,60]. 

A rigorous analysis of these models can deepen understanding of 
the translation process, assist in integrating the vast amount of 
empirical findings, lead to efficient algorithms for optimizing gene 
translation for various biotechnological goals, and help to improve 
the fidelity and predi( ti\'(^ ability of the models. In the near future, 
this will enable building syntectic biological devices that are based 
on re-engineering biological mechanisms and specificciUy gene 
expression. 

A newly developed technique, called ribosome profiling [61,62], 
provides indications on the occupancy of codons by ribosomes 
along the mRNA molecules in vivo. This breakthrough has led to 
a renewed interest in computational models for translation (see 
e.g., [63-65]). 

Reuveni et al. [66] considered a deterministic model for 
translation called the ribosome flow model (RFM). This model is a 
deterministic approximation of an important model from statistical 
physics, called the asymmetric simple exclusion process (ASEP), that is 
the standard mathematical model for ribosome flow. ASEP has 
also been used to model and analyze many other systems and 
processes, including traffic flow, molecular motors, surface growth, 
the movement of ants, and more [67]. 

In this paper, we study the dynamical behavior of the RFM 
under the assumption that some or all of its parameters vary 
periodically, with a t:ommon period T. This models periodically 
time-varying initiation and/or transition rates along the mRNA. 
We refer to this model as the periodic: ribosome flow model (PRFM). 

Main results and their implications. Our main result 
shows that the PRFM entrains to a periodic excitation. In other 
words, the PRFM admits a unique periodic solution with period 
T, and all the state-variables converge to this solution. This means 
that all the ribosome densities converge to a periodic pattern with 
period T and, in particular, the protein synthesis rate converges to 
a periodic pattern. To the best of our knowledge this is the first 
proof of entrainment in a non-trivial mathematical model for 
translation. 

Our results suggest that entrainment takes place in particular in 
the case where the codon decoding rates (called transition rates in 
the RFM) are constant, and the initiation rate is T-periodic. 
Similarity, entrainment takes place if the initiation rate is constant 
and some of the transition rates are T-periodic. From a biophysical 
perspective, this suggests that periodic oscillations of the transla- 
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tion rate (and thus protein abundance) can be induced in various 
ways including: 1) oscillations of factors related to the initiation 
step such as the mRNA levels of genes, the abundance of 
ribosomes, and the abundance of initiation factors; and 2) 
oscillations of factors related to the elongation step such as the 
abundance of c-longation factors and tRNA genes. Specifically, 
oscillations in the abundance of a single tRNA gene is enough to 
induce oscillations in the translation rate and protein abundance. 

These results have several implications. First, they support the 
conjecture that cell-cycle dependent dynamics of proteins may be 
obtained by entrainment in the translation process. Moreover, the 
biological mechanism can generate a periodic production rate 
relatively easy; it is enough to vary just one tRNA abundance in a 
periodic manner. However, in the PRFM entrainment takes place 
whenever the initiation/transition rates vary periodically (with a 
common period) regardless of their amplitude. This suggests that the 
bottleneck argument in [12] is not necessarily needed. 

Second, in the context of synthetic biology our results may lead 
to new mechanisms for generating various syntectic devices at the 
translation level. Several recent studies considered the design of 
synthetic biological oscillators, mostly based on manipulating 
aspects related to transcription (see, e.g. [1,8,68-71]). The authors 
of [13] raise the question of why would cells regulate translation 
using codon usage and changes in tRNA modification status. They 
hypothesize that a rapid change in the abundance of tRNA 
modifications may allow cells to quickly reset the translation speed 
of existing transcripts, and thus respond quickly to stress or other 
changes in environmental conditions. If this is indeed so then 
developing synthetic biology devices based on entrainment at the 
translation level may have unique advantages. 

Mathematical tools. In order to make this paper more 
accessible, we now briefly explain the main mathematical tools 
that are us(;d in analysis. 

Proving entrainment in non-linear dynamical systems is non- 
trivial. One standard approach is based on contraction theory 
[3,72,73]. A dynamical system is called contracting if the distance 
between trajectories emanating from any two initial conditions 
quickly decreases with time (more precisely, it decreases at an 
exponential rate). This means that the information about the 
initial condition is "quickly forgotten". 

Consider a system that is periodically excited with a period T. 
Assuming that the trajectories remain bounded, it is possible to 
show that the system admits a periodic solution with period T. 
Consider two trajectories, one emanating from an initial condition 
on this periodic solution, and the second from some arbitrary 
initial condition. If the system is also contracting then these 
trajectories must converge to one another, so all trajectories 
converge to the periodic solution. This proves entrainment. 

The proof of entrainment in the PRFM is based on these ideas. 
However, some additional analysis is needed, as the RFM is on the 
"verge of contraction", yet it is not contracting on its entire state 
space. 

The remainder of this paper is organized as follows. The next 
section briefly reviews the ASEP and RFM. The main results 

about entrainment in the PRFM are described in the Results 
section. The proofs are detailed in the Methods section. The 
Discussion section provides a summary, and describes possible 
directions for further research. 

Preliminaries: From ASEP to tlie RFIVI 

An important computational model for translation elongation is 
the Asymmetric Simple Exclusion Process (ASEP) [25,26,74]. In this 
stochastic model particles hop, according to some probability 



function, between consecutive sites on a ID lattice. Each site can 
be either occupied by a particle or not. Hops may take place only 
to a target site that is not already occupied by another particle 
(hence the term simple exclusion). The term asymmetric implies that 
there is some preferred direction of movement along the lattice. If 
motion is allowed only in one direction then ASEP is sometimes 
called the totally asymmetric simple exclusion process (TASEP). ASEP 
(and its many variants) is regarded as a paradigmatic model for 
non-equilibrium statistical mechanics and has been used to model 
and analyze various biological systems and processes, including 
intracellular transport, molecular motors, pedestrian dynamics and 
of course gene expression [75-80] . 

In TASEP models for translation each site in the lattice 
corresponds to a codon, the hopping particles are ribosomes, and 
their footprint includes £ sites (in the case of translation, the 
footprint of a ribosome is usually 9<£< 12 codons corresponding 
to 28 — 35 nt [61]). For example, a new particle (ribosome) can 
enter the lattice only if all the first £ sites are all empty. Initiation 
time as well as the time a ribosome spends translating each codon 
are random variables [e.g., with an exponential distribution), and 
are codon-dependent. Analysis of TASEP is based on determining 
the probability of steady-state configurations using matrix products 
(see the excellent survey paper [81]). 

Reuveni et at. [66] recently considered a deterministic model for 
translation called the Ribosome Flow Model (RFM). The RFM is a 
finite-dimensional mean-field approximation of TASEP (see, e.g., 
[81], p. R345 and [82], p. 1919). The RFM includes n state- 
variables connected by a set of n non-linear first-order differential 
equations: 

x\ =X{\ —x\) — k\X\{\ —X2) 
i;2 = AiXi(l -X2)— 22^2(1—^3) 
±3 = ^2x2(1 — x-}) — 2.3x3(1 —X4) 

(1) 

Xfj — i = An— 2Xk— 2(1 ^n—l) 2ff—iXfi—i(\ Xjj) 
Xfi = Afj — i Xfi — 1 ( 1 Xj^ 2.nXfi . 

The positive parameters X . . . ,X„] are called the initiation rate 
[transition rates]. The state-variable x,- : ->[0,1] describes the 
"level of occupancy" of ribosomes at site i at time t, where 
Xi(t) = 1 [Xi{t) = 0] corresponds to the site being completely fuU 

[empty]. 

To explain this model, consider for example the equation for Xi , 
i.e. the change of "level of occupancy" at site I. The term —xi) 
models the fact that ribosomes reach the first site with initiation 
rate 1, but their effective binding rate depends on how occupied 
site 1 is. In particular, if Xi(f)= 1, i.e. the site is completely full, the 
effective binding rate is zero. The term — 11X1(1 —X2) describes 
the rate of transition of ribosomes from site I to the consecutive 
site 2. 

Note that in ASEP each site may include either zero or one 
particles. In the RFM, the X/s correspond to averaged occupancy 
levels and therefore x,(?) takes values in the closed interval [0,1]. 

We refer to 

R(t) : =X„x„itl 

that is, the exit rate of ribosomes from the last site, as the translation 
rate at time t. 
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Since the state-variables correspond to normalized occupancy 
levels, the initial condition x(0) is always in the closed unit cube: 

C: ={xeK":Xi6[0,l],/=l,...,n}. 

The simulation results in [66] show that TASEP and its mean-field 
approximation (the RFM) yield similar predictions of translation 
rates. For example, the correlation between their predictions over 
the set of cutlogcnous genes of S. cerevisine is 0.96. Important 
features of translation elongation that are captured in TASEP, for 
example, the sequential order of the codons, translation eflSciency, 
the interaction between ribosomes and their jamming, the 
initiation and elongation rates, are also encapsulated in the RFM. 

We now briefly summarize some known results on the 
dynamical behavior of the RFM. Let 

Int(C) : ={xe]g": x,€(0,l),/= 1, . . . ,«} denote the interior of C, 
and let x{t,d) denote the solution of the RFM at time t for the 
initial condition x(0) = aeC. It has been noted in [83] that the 
RFM is a (tridiagonal) monotone dynamical system [84]. Combining 
this with the fact that C is an invariant set of the dynamics and a 
theorem of SmUlie [85] yields the following result. 

Theorem 1 \83\ The RFM admits a unique equilibrium point 
eelnt(C), and lim,_,oo x(t,a) = efor all aeC. 

This means that there exists a unique steady-state profile of 
ribosome distributions (and thus a unique translation rate). The 
trajectory starting from any initial distribution will com erge to this 
steady-state profile. Changing the values of the positive parameters 
will not change this qualitative picture, only the steady-state 
profile. 

Let I 'I] : denote that Li vector norm, i.e. 

kli = kl|+ ■•• +|"k|. It has also been shown in [83] that the 
RFM is non-expanding with respect to the Li norm, that is, 

\xit,a)-xit,b)\i<\a-b\i (2) 

for all t>0 and all a,beC. This means that the Li distance 
between two ribosome distribution profiles can never increase. It is 
worth noting that both Theorem 1 and (2) follow immediately 
from the more general results in this paper. 

In some cases the transition rate along genes is constant [86] , so 
the translation eHiciencies of all the codons are identical. This 
happens, for example, when the rate limiting factor is the 
concentration of elongation factors and not the local features of 
the coding sequence, such as tRNA molecules or when there is a 
balance between the codon frequencies and tRNA levels [87]. In 
the context of the RFM, this can be modeled by considering the 
special case where 

= ^2 = • ■ • = = : Ac, 

that is, the transition rates A, arc all e(|ual, and Ic denotes their 
common value. Since this Homogeneous Ribosome Flow Model 
(HRFM) includes only two parameters, X and Ac, the analysis is 
considerably simplified. Ref. [88] analyzes the qualitative and 
quantitative dependence of e on the parameters /l,Ac in the 
HRFM. Ref [89] studied the HRFM when w^oo, i.e. when the 
length of the mRNA chain goes to infinity. In this case, it is 
possible to obtain closed-form expressions for the equilibrium 
point using the theory of infinite 1 -periodic continued fractions. 

In eukaryotes the translation rate can affect the initiation rate 
via recycling of ribosomes. To model this, Ref [90] has considered 
the RFM as a control system. Here 1 is replaced by a function u(t) 
(the input), and an output y{t) = R{t) is added. It has been shown 



that this is a monotone control system, as defined in [91]. Also, analysis 
of the closed-loop system, obtained by closing the loop from ytou 
with positive linear feedback, has shown that several nice 
properties of the RFM hold also for the closed-loop system. In 
particular, there exists a unique globally asymptotically stable 
equilibrium point e. For the special case of equal A,s, closed-form 
expressions relating the closed-loop system parameters and e have 
been obtained. 

Results 

A function / : ]R^]R is called T-periodic if 

m=f(t+T) (3) 

for all t. For example, sin (?) is Ink periodic, for all integers k. In 
this paper, we consider the behavior of the RFM (1) under the 
following assumptions (see Figure 1): 

1. The initiation rate X(t) and transition rates h(t), /=!,...,«, 
are continuous, strictly positive and uniformly bounded 
functions of time, i.e., there exist ^i,i52£lR such that 

0<di< Ai(0,A(0 < Si, for all ? > 0. (4) 

2. There exists a minimal T>0 such that all these functions are 
T-periodic. 

We refer to this case as the periodic ribosome flow model (PRFM). 

Remark 1 Note that the PRFM includes in particular the case 
where some of these rates are constant, as a constant function is T- 
periodic for every T . However, item 2) above implies that the case 
where all the rates are constant is ruled out, as then the minimal T 
is zero. Indeed, this case is just the RFM. 

The next example illustrates the dynamical behavior of the 
PRFM. 

Example 1 Consider the PRFM with « = 3, 
A(0 = 2-l- sin(27c0, 



;i2(0 = 3-l- sin(27c?-l-l/2), 



A3 (0 = 4 -h 2 cos (47cf -h 1 /8), 



and initial condition x(0) = [ 1/2, 1/2, 1/2 ]'. Note that all the 
rates here are periodic, with minimal common period T=\. 
Figure 2 depicts Xi(t), i = 1,2,3, as a function of t. It may be seen 
that each state-variable converges to a periodic signal with period 
T=\. 

The next example considers the case where only the initiation 
rate oscillates, and the constant transition rates are the rate- 
limiting factors. 
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Figure 1 . Upper part: the elongation rates of codons and the initiation rate are T-periodic, for example, due to signals related to the cell-cycle. Lower 
part: in the RFM, this is modeled by T-periodic rates X{t) and yielding the PRFM. Our main result shows that consequently the translation rate 
and ribosomal densities (the a-,s) converge to a unique T-periodic solution. 
doi:1 0.1 371 /journal.pone.0096039.g001 



Example 2 Consider the PRFM with n = 3, 



A(r) = 2+ sin(77:/), 



0.8 




0.4 -: V 

\ -, 




0 0.5 1 1.5 2 2.5 3 3.5 4 



t 

Figure 2. State-variables x^lf) [solid line]; xzifl [dashed]; and 
xjift [dotted] (y-axis) as a function of time (x-axis) in Example 1. 

All state-variable converge to a periodic signal with period T=l. 
doi:10.1371/journal.pone.0096039.g002 



and initial condition x(0) = [ 1/2, 1/2, 1/2]'. Note that here 
the transition rates are constant and relatively small, whereas the 
initiation rate is larger and periodic with period T = 2. Figure 3 
depicts Xi(t), i = 1,2,3, as a function of?. It may be seen that all the 
x,s converge to a periodic signal with period T = 2. Note that the 
rate-limiting transition rates considerably attenuate the oscillations 
amplitude as they propagate through the mRNA chain. Here the 
small transition rates form the "botdeneck" in the process. This 
example demonstrates that entrainment takes place even when the 
rate-limiting factor is constant, as long as there exists at least one 
other factor that is periodic with T>0. 

Our main result shows that the state-variables in the PRFM 
always entrain to a unique periodic solution. Let x{t,to,a) denote 
the solution of the PRFM at time t for the initial condition 
x(to) = a. 

Theorem 2 The PRFM admits a unique periodic solution 
y ■ M+ "'■IritCC), with period T, and 
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0.9 




0.4 - •-. 

0.3 1 1 < 

0 5 10 15 

t 

Figure 3. State-variables x^if) [solid line]; X2{fl [dashed]; and 
xjift [dotted] (y-axis) as a function of time (x-axis) in Example 2. 

The initiation rate is periodic with period r = 2, while the transition 
rates are constant and relatively small. All state-variable converge to a 
periodic signal with period T = 2, but the amplitude of the oscillations 
is considerably attenuated as it passes through the mRNA chain. 
doi:1 0.1 371 /journal.pone.0096039.g003 



lim |x(/,0,fl)-y(?)li =0, for all aeC. 



In other words, every trajectory converges to the unique periodic 
trajectory y. In particular, the translation rate R{t) = A„(t)Xii{t) 
converges to the T-periodic function A„{t)y^{t). 

By Remark 1, Theorem 2 holds in particular in the case where 
the transition rates Xj, i=l, . . . ,n, are constant, and the initiation 
rate is T'-periodic (but not constant). Similarity, Theorem 2 also 
holds if the initiation rate is constant and some of the transition 
rates, A,, are 7"-periodic (but not constant). 

The stochastic nature of reaction events induces random noise 
in biochemical networks (see, for example, [92,93]). As noted in 
[94], this becomes particularly important when there are few 
molecules in the system, as is often the case in a cell. It is natural to 
consider whether entrainment in the PRFM takes place also in the 
presence of noise. Our simulations suggest that this is indeed the 
case. 

Example 3 Consider the PRFM with « = 3, 
Mt) = 2+ sin(27iO + vfi(0. 



;.i(0=i+tt'2(0, 



X2{t) = 3+ sin(27i;-|-l/2) + u'3(/), 



i3(/) = 4 + 2cos (47r/+l/8) + H'4(r), 

and initial condition x(0) = [l/2, 1/2, 1/2]'. The U',(?)'s are 
drawn as independent random values from the uniform distribu- 



0.8 




t 



Figure 4. State-variables xi((t [solid line]; X2{t) [dashed]; and 
X3(ft [dotted] (y-axis) as a function of time (x-axis) in Example 3. 

The initiation and transition rates are periodic with a common period 
r= 1, but with added random noise. It may be seen that each state- 
variable converges to a periodic signal with period T=l, but with 
added noise. 

doi:10.1371/journal.pone.0096039.g004 

tion on the interval ( — 1,1) for all t. Note that this implies that all 
the rates remain positive for all /. Figure 4 depicts /= 1,2,3, 

as a function of t. It may be seen that each state-variable still 
converges to a periodic signal with period T=l, but with 
perturbations induced by the noise. 

Further study of entrainment in the PRFM in the presence of 
random perturbations is beyond the scope of this paper. However, 
we note that there exist theoretical results on contraction in the 
presence of random noise; see e.g. [95,96]. 

As mentioned above, the RFM is a mean-field approximation of 
TASEP. Thus, our results suggest a natural question, namely, does 
TASEP entrain? 

Example 4 Consider a TASEP model with sites. When a 
particle in site i cannot hop because site i+l contains a particle, 
the next hopping time is determined by 

where t is the current time, and p(li) is a random variable drawn 
from the exponential distribution with mean parameter Recall 
that the exponential probability distribution function with mean m 
is given by 

fix): = —exp( — x/m), forallx>0. 
m 

At the next hopping time, this particle hops unless the next site is 
fuU again, in which case a new hopping time is drawn. 

We ran a simulation of this process with A^= 10, and final time 
Tf = 200000. In the simulation, performed using MATLAB, time 
was discretized using a time step of 0.01. (In particular, the next 
hopping times are always rounded to a value k * 0.01, where k is 
an integer.) Initially, all sites are empty. The rates are 
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^2(0 = 3+ cos(^ + 



2nt 
20000 



13(0 = 2+ sin(l + J^), 

and 1,(0 = 2 for all i^{2,3}- Note that all these rates are periodic 
with a common minimal period 



r = 20000. 



(5) 



Fig. 5 depicts the results. The time range [0,7/] is divided into 
segments of length 1000 (so that there are r//1000 = 200 
segments), and the 0/1 occupancy at each time step is averaged 
on each segment. For example, the value depicted at time segment 
1 is the averaged occupancy in the time interval [0,999]. The 
averaged occupancy on each segment is shown in Fig. 5 for sites 
2,3,7, and 9. It may be seen that aU the averaged occupancies 
entrain to the periodic excitation. In particular, they are periodic 
(up to the noise induced by the stochastic process) with a period of 
20 segments, corresponding to a time period of 20 * 1 000 which is 
equal to T in (5). 

Our simulations do suggest that some form of entrainment also 
takes place in TASEP. Of course, one must first rigorously define 
what entrainment means in a stochastic model such as TASEP. 

Discussion 

Many biological and physiological processes are periodic (see, 
e.g., [97]), indicating periodicity in their corresponding gene 
expressions. For example, the cell-cycle is a periodic series of 
events that allows cells to replicate. 

Two recent studies suggest that the protein levels of ceU-cycle 
related genes are regulated by periodically varying tRNA levels 
[12,13]. In other words, the translation-elongation mechanism 
entrains to these periodic oscillations. 

To examine the plausibihty of this idea, one must first consider 
the time constants involved; specifically, the mRNA life time and 
translation time should be longer than the cell-cycle period. For 
concreteness, consider the case of .S*. ceremsiae. The cell-cycle period 
in S. cerevisiae is less than 87 minutes [98] (ceU-cycle period can be 
much shorter in eukaryotes; for example, it was reported that the 
duration of ceU-cycle in early embryo of the fruit fly D. melamgaster 
is only 8 minutes [99]). In S. cerevisiae there are hundreds of genes 
with mRNA half-life larger than 100 minutes (see, for example, 
[100]). The translation rate in S. cerevisiae was estimated to be 
higher than 0.956 codons per second (the slowest codon is CUU) 
[101] with average rate over all codons of 10 codons per second (in 
mouse the average codon translation rate was estimated to be 
around 5 codons per second [86]). In practice, this rate can be 
much slower due to strong folding of the mRNA and interaction of 
the translated amino-acids with the exit channel of the ribosome 
[27], In S. cerevisiae the ORE length range is between 51 and 
14,733 nucleotides; for example, the longest gene in yeast is 
MDN1/YLR106C which includes 4,91 1 amino acides (see http:// 
www.yeastgenome.org). This is a huge dynein-related AAA-type 
ATPase (midasin) which forms extended pre-60S particle with the 
Rixl complex (Rixlp-Ipilp-IpiSp) and acts in removal of 
ribosomal biogenesis factors at successive steps of pre-60S 
assembly and export from nucleus. This gene corresponds to an 
upper bound on the translation time of a gene that is larger than 8 1 




200 



# time segment 



Figure 5. Simulation of TASEP with periodically time-varying 
rates. The plot shows the averaged occupancy over time segments as a 
function of the time segment. 
doi:10.1371/journal.pone.0096039.g005 



minutes (assuming a lower bound on translation rate of 1 codons 
per second; which may be lower in practice); an estimated 
translated time of this protein based on mean codon translation 
time is 8.2 minutes. In mammals the mean codon decoding time is 
5 codons per second and the longest human protein (Titin — TTN) 
which has 33,000 amino acids, corresponding to estimated 
translation time of 1 1 0 minutes. This suggests that periodically 
varying tRNA levels may indeed induce periodic expression levels 
of cell-cycle related proteins. In addition, assuming that time to 
steady sate is related to the translation time (at least one ribosome 
should finish the translation), an estimated lower bound on the 
oscillations is 8/1 10 a; 1/13.7 the time to reach steady state. 

To rigorously analyze entrainment at the translation level, we 
considered the REM under the assumption of periodic initiation 
rate !(/) and/or periodic transition rates l/(?) with a common 
minimal period T>Q. Our main result is that all the ribosome 
densities converge to a unique periodic solution with period T. 
This implies in particular that the protein translation rate 
converges to a unique periodic function. The PRFM is thus the 
first computational tool providing an explanation of how 
periodicity can be passed from the translation to the protein level 
via the codon usage bias, i.e., the differences in the frequency of 
occurrence of different codons in the coding sequence. Specifical- 
ly, according to the PRFM the distribution of codons in different 
open reading frames (OREs) can affect their oscillations. For 
example, if the levels of certain tRNA species oscillate this wiU 
affect only genes with codons that are recognized by these tRNAs. 

These results support the conjecture that oscillations of the 
tRNA levels and/ or initiation factors, with a common period T, 
induce periodicity in the protein levels of ceU-cycle genes. The 
assumption of variations with a common period may seem 
unjustified, but analysis of the PRFM shows that it is enough to 
oscillate a single tRNA level, or just the initiation rate to obtain an 
oscillatory behavior, with the same period, in protein synthesis. 
Furthermore, genes that are part of a certain pathway and/or 
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function usually have common regulators (see, for example, [102]), 
and this may lead to periodic oscillations with a common period. It 
is important to emphasize that we do not claim here that the 
oscillations amplitude is necessarily "large". There may be cases 
where the amplitude of the oscillations in some x, or R may be 
small and in this case the signal may "look" constant; for example, 
in the case that the oscillations are not in the botdeneck (in terms 
of translation rate) of the gene they may have a smaller effect on 
the translation rate. 

Oscillations and entrainment also play an important role in 
synthetic biology- [1,8,68—71]. Ind(;ed, a major challenge in this 
field is scaling up to larger and more complex biological systems. 
One possible approach is to design networks based on an 
interconnection of several biological elements (modules) that 
synchronize to a single central clock. Entrainment is needed to 
achieve this. In this context, the period T of the oscillator may 
perhaps be controlled and, in particular, made much shorter than 
the cell-cycle period. The analytic results on the PRFM may thus 
lead to new synthetic devices that produce periodically-varying 
protein levels based on oscillations in tRNA levels and/or 
initiation factors. Indeed, the analysis of the PRFM suggests that 
there are many different possible ways for generating such a 
periodic dynamics. 

It is important to remember that oscillations in various factors, 
and not only the concentrations of tRNA molecules, can affect the 
periodic dynamics of the translation process. Among others, the 
oscillations in the concentrations of Aminoacyl tRNA synthetase, 
ribosomes (via changes in concentrations of ribosomal RNA genes 
and/or ribosomal proteins), elongation and initiation factors, 
mRNAa levels, and free amino acids may trigger an oscillatory 
behavior of the translation process. 

One experimental approach for validating our theoretical 
results, and for further research of oscillations in translation is by 
in-vitro single-molecule fluorescence experiments (see, for example 
[103-107]). Such an experiment should encompass the different 
components of the translation machinery (tRNA molecules, 
ribosomes, elongation and initiation factors). Specifically, it will 
enable monitoring oscillations at the level of single ribosomal 
movements. 

There are other types of large scale in-\ itro c'xperiments (see for 
example, [108-111]) that may be relevant. Here, crude extracts 
containing all the macromolecular components (70S or 808 
ribosomes, tRNAs, aminoacyl-tRNA synthetases, initiation, elon- 
gation and termination factors, etc.) required for translation of 
RNA are prepared. To ensure efficient translation, each an extract 
must be supplemented with amino acids, energ)- sources (ATP, 
GTP), energy regenerating systems (creatine phosphate and 
creatine phosphokinase for eukaryotic systems, and phosphoenol 
pyruvate and pyruvate kinase for the E. coli lysate), and other co- 
factors (Mg2+, K+, etc.). 

In-vitro experiments should enable to produce oscillations in 
different molecules related to the translation process (for example, 
the tRNA levels; similarly to the mathematical example described 
above) and measuring the effect of these oscillations on the 
ribosomal translation pattern. 

The mathematical analysis performed here leads to several 
computational questions that deserve further study. First, our 
results provide little information on the periodic solution y. In 
particular, important questions are how does the amplitude and 
other properties of y depend on the initiation and translation rates, 
and what is the convergence rate of the solutions to y. The analysis 
suggests that there are several possible ways to induce periodicity 
in the translation rate (e.g., via oscillating tRNA molecules, mRNA 
molecules, initiation factors, elongation factors, ribosomal RNA, 



etc), and it would be interesting to analyze how the periodic 
solution y is affected by the different possible ways of inducing 
periodicity. 

Finally, and more generally, TASEP and its variants have been 
used to model and analyze a large number of biological and 
artificial systems including biomolecular motors [112,113], the 
collective motion of ants [114], traffic flow [11, 5- 117], ad hoc 
communication networks [118], and surface growth [119]. Many 
of these systems may be affected by periodic signals. For example, 
traffic flow is often controlled by periodically-varying traffic lights. 
It may be of interest to model and study such systems using the 
PRFM. 

Methods 

In order to prove our main results, we first detail several known 
results that will be used later on. 

Preliminaries 

Consider the system 

x=fit,x), (6) 

evolving on a convex set Cc]^". Let x(t,to,Xo) denote the 
soluticm of (6) at time t>tQ with x{to) = Xo (for the sake of 
simplicity, we assume from here on that x{t,tQ,XQ) exists and is 
unique for all t>to>0). 

Recall that (6) is said to be contraeting [72] on C with respect to a 
norm |-| : if there exists OO such that 

|xfe,fi,a)-x(f2,fi,^)|< exp(-(f2-fi)c)|a-6| (7) 

for all t2>ti>0 and all a,beC. In other words, trajectories 
contract to one another at an exponential rate. 

A standard approach for proving contraction is based on 

analyzing the Jacobian matrix J(t,x) : = — f{t,x). Recall that a 

vector norm |-| : induces a matrix norm \\-\\ : M"~*M+ 

defined by 

Pll : =max|^x|, 
1*1 = 1 

and a matrix measure n : defined by 

H(A): =lim-(||/-h€^||-l). 

The next result, known as CoppeFs inequalitv (see, e.g. [120,121]), 
provides a bound on the solution of a linear time-varying system in 
terms of the induced matrix measure. 
Theorem 3 Consider the differential equation 

m = M(t)y(t) (8) 



where >' : ]K+ ~>M"j ^"-^ the matrix M : ->]R"^" " defined and 
continuous for all t>to>0. Suppose that |-| : ]K"— vector norm 

and let n : ^ " denote the induced matrix measure. Then every 
solution of (8) satisfies 
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\y{t)\ < exp 



li{M{s))ds ) b(ro) I , for all t>to. 
to J 



We now give an informal explanation of how this can be used to 
prove contraction. Suppose that there exists c>0 such that 



H{J{t,x)) <—c 



(9) 



for all t>0 and all xeC, where C^J^" is a convex set. The 
distance dx between two trajectories of (6) emanating from 

. d ^ 

infinitesimaUy close initial conditions satisfies ^ ( 3x) = J{t,x) Sx, 

so combining (9) and Coppel's inequality suggests that (7) holds for 
all a,beC. Furthermore, it can be shown that contraction implies 
entrainment to a periodic excitation, see e.g. [3,72] for rigorous 
statements and proofs. 

Consider applying these ideas to the PRFM. We can write the 
PRFM as x=f{t,x), where / is T-periodic, \.e.f(t+T,x)=f{t,x) 
for all X and t. A calculation shows that the Jacobian of / is 



J{t,x) = L(t,x)-D{t), 
where L is shown below and 

D: = diag(A,0, 0,A«). 



— X\{\—X2) XlXl 
Xl(\—X2) — XlXl — X2(\ — Xl) 
0 l20--Xi) 



(10) 



L = 



— Xn-lXn-1 — 

Xn-\(y ~Xn) 
X„-l(l-X„) 



0 

0 
0 

Xn—lXn—l 
— Xn — lXn-l 



Recall that a matrix is said to be a Metzler matrix if all its o£F- 
diagonal entries are non-negative. Note that L(t,x) is Metzler, 
tridiagonal, and with zero sum columns for all ?>0 and all xeC. 

It is well-known ([122], Chapter 3) that the induced matrix 
measure corresponding to the L\ vector norm is 



where 



Hi{A)= max{ci(^), . . . ,c„{A)}, 



l<i<n 



(11) 



(12) 



i.e., the sum of the entries in column /, with non diagonal elements 
replaced by their absolute values. Of course, if ^4 is Metzler then 



one can take Ay instead of \Aij\ 



(12). 



Calculating /ij {J{t,x)) for the PRFM shows that when n = 2, 

jll{J{t,X)) = msx{ — k(t),— A2{t)} 



so in this case we have contraction with respect to the L\ norm. 
However, for n>3, c,(/(;,x)) = 0 for i = 2,3, . . . ,n— 1, and 
Hx{J{t,x)) = Q. Intuitively, this means that the PRFM is on the 
"verge" of contraction with respect to the Li norm, but this is not 
enough to prove entrainment. Furthermore, when X\=Q and 
X3 = 1 all the entries in the second column of / are zero, and this 
implies that the PRFM is not a contraction on C with respect to any 
norm (as a necessary condition for contraction is that J{t,x) is a 
Hurwitz matrix for all f>0 and all xeC). 

Proof of Main Result 

For two vectors a,be^, we write a<Z) if a, <Z),- for /= 1, . . . ,n. 
Let InSjR" denote the vector with all entries equal to one. For 
Ce[0,l/2], defme 



= {xeC:a«<x<(l-01«}- 



Note that Qq = C, and that is a strict subcube of C for ail 
Ce(0,l/2]. The next result shows that the trajectories of the RFM 
always enter such a strict subcube, and then remain in it. 

Proposition 1 Consider the PRFM. Fix arbitrary t\ > 0, aeC, and 
T>0. There exists C = C(T)e(0,l/2), tuith C(t)->0 as t->0, such that 



x(?,ii,a)eQj, forallt>t\ + x. 



(13) 



Proof. See the section Additional Proofs. 

Recall that the PRFM is not a contraction on C with respect to 
any norm. The next result shows that contraction does hold on 
any strict subcube of C. The proof, given in the Additional Proofs 
section, uses a suitable diagonal scaling of the L\ norm. 

Proposition 2 Fix an arbitary fe(0,l/2). Then the PRFM is 
contracting on Q.^. 

Note that all the proofs up to this point (given in the section 
Additional Proofs) do not rely on the assumption that the rates are 
periodic functions. We can now prove our main result. 

Proof of Theorem 2. Recall that the excitation is periodic with 
period T. Let m(a) : = x{T,0,a), i.e., m maps the initial condition 
x(0) = « to the solution of the PRFM at time T. Then m is 
continuous and maps C to C, so by the Brouwer fixed point 
theorem (see, e.g. [123]) there exists t^eC such that w(0 = C, i.e. 
x(7',0,C) = C- This implies that the PRFM admits a periodic 
solution y with period T. It follows from Proposition 1 that 
x(?,0,Q6lnt(C) for all ? > 0. We already know that every trajectory 
enters some strict subcube of C, that this subcube is an invariant 
set, and that in this subcube contraction holds. Thus, every 
trajectory converges to the periodic solution emanating from 
x(0) = C- In particular, there cannot be two distinct periodic 
solutions. This completes the proof of Theorem 2. 

Note that the reasoning above does not rule out the possibility 
that 



y(0=y(0), 



(14) 



for all / > 0, i.e. that the periodic trajectory is just a fixed point. To 
study this case, assume that a fixed point esC indeed exists. Then 
(1) yields 

A(0(l-ei) = Ai(0ei(l-e2) 
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= /l2(0e2(l-e3) 



= A„_i(Oe„-i(l-e„) 



= A„it)e„, 

for all t>0. Thus, (14) is possible only in the rather special case 
where all the rates are equal, up to multiplication by a non- 
negative scalar. 

Additional Proofs 

This section includes the proofs of several results stated above. 

We begin by .stating and proving two auxiliary results. The next 
subsection describes a result that wiU be used to prove Proposition 
1. 

Repelling Boundaries and Persistence 

Lemma 1 Consider a time-varying system 

x=fit,x) (15) 

evolving on a subset ofX: =I\ xlix ... x /„ c J^*^ , where each Ij is 
an interval of the form [0,A], A>0, or [0,oo). Suppose that the time- 
dependent vector field f =\f\, . . . f„ \ has the folhwing boundary-repelling 
property: 

(BR) For each (5 > 0 and each sufjicientlj small A > 0, there exists 
K = K{5,/^>Q such that, for each k=l, . . . ,n and each ?>0, the 
condition 

Xk<A, and Xi>5, for every \<i<k—\ (16) 

(for k=\, the condition is simply x\ < A) 
implies that 

fk(t,x)>K, forallt>0. (17) 



Then given any t>0 there exists £=£(t)>0, withe{z)^0 t->0, such 
that, for every solution x{t), f>0, it holds that 

Xi{t)>e for all ...,«} and all t>%. 



In other words, the contJusion is that after an arbitrarily short time 
every Xi{t) is separated away from zero. 

Proof of Lemma 1. Pick any T>0. Let T„ : =xln. We proceed by 
induction: for each /re{ 1, . . . ,«} we will define an Ek and show that 
for every solution x(f), Xi{()>ek for every t>kx„ and every 
ie{l, . . . ,/:}. Then £=£„ gives the result. Pick any frxed A>0 
small enough that (BR) holds, and let x{t) be any given solution. 
From here on, we write just T instead of T„. 



Consider first the case fc = 1. Let K = K(^t^,l)) (for any arbitrary 
(5 >0) and define £i : = min{A,.SrT}. Let ?oe[0,T] be such that 
xi(?o)^£i- Such a ?o exists, since by property (BR), x\{t)<e\ <A 
for all ?e[0,T] would imply that X\{t~)=f\{t,x{t~))>K>Q for all 
?6[0,t], which in turn implies X\{x)'>X\(^)^Kx>Kx>e\, con- 
tradicting xi(t)<£i. We claim that also Xi(;)>£i for every ?o. 
Indeed, suppose otherwise. Then, there is some t\ > to such that 
( : =xi(ri)<£i. Let 

a: = min{t>to\xi{t)<i}>to. 

As Xi((7)<f<£i < A, property (BR) says that 
Xi{(j)=fi{a,x(a))>K, so it follows that x\(t)>0 on an interval 
[f7 — t,(t], for some t>0. But then X]{a — T)<Xi{(j)<£, which 
contradicts the minimality of (7. Thus Xi(f)>£i for all t>to, and 
in particular for all t>z. 

Now by induction, consider l<k<n, and suppose that 
Xi{i)>£k for all t>kT, and every ie{l, . . . ,k}. We must define 
Sk+i so that Xi(t)>ek+\ for all />(/r+l)T and every 
(e{l,...,A:-|-l}. Let K = K(ek,A) in (BR), and define 
£k+i ■■ = rmn{Sk,A,KT:}. 

Let toe[kx,{k+l)T:] be such that Xk+i(to)>ek+\. Such a to 
exists, since by property (BR), (using that Xi{t)'>£k= : d for all 
i<k and t>kx) Xk+\(t)<£k+i <A for all te[kx,(k-\-l)x] would 
imply that Xk+i(t)=fk+iit,x(0)^K>0 for all te[kx,{k-\-l)x], 
which in turn implies 

Xk+i((k-\- l)x)>Xk+i(kx) + Kx>Kx>£k+i, contradicting 
Xk+i{(k+l)x)<ek+i. 

We claim that also Xk+i(t)>£k+\ for ever)- t>to. Indeed, 
suppose otherwise. Then, there is some tk+i>to such that 
i ■ =Xk+iitk+i)<£k+l- Let 

a: = mm{t>to\xk+\(t)<^}>to■ 
As xt+i(iT)<^<£t+i <A, and also Xi{(T)>£k (because Xi(t)>£k 
for all t > kx, and (T>to> kx), we may apply property (BR), which 
says that Xk + ]((r)=fk + ](cy,x(ff))>K, so it follows that Xk+\(t)>0 
on an interval [ff — t,(t], for some T>0. But then 
Xk+\{''^i^)<Xk+\i(>)<i, which contradicts the minimality of 
(7. Thus Xk+i(t)>ek+i for all t>to, and in particular for all t>x. 
By the definition ot Ek+l and the induction hypothesis, we also 
have Xi{t)>Sk+l for all t>(k+l)x and every fe{l, . . . ,k}. 

We can now prove Proposition 1 . We begin by showing that the 
PRFM satisfies property (BR) in Lemma 1 on X = C. Fix an 
arbitrary ^ >0. Consider the case k=\. If Xi <A then 

/l=/l(l-Xi)-^lXi(l— X2) 

>Ku 

where .fiTi : = S\(l — A)— SiA. Now pick 1 <fc<n. If Xi<A and 
Xi>d for 1</<A:— 1 then 

fk = h-lXk-iil-Xk)-hxkil-Xk+i) 
>K, 

where 
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1. Hi ^L)= max{ci,C2, . . - ,c„} where 



K : = did{l-A)- diA. 
Finally, if x„<A and Xi>d for l</<n— 1 then 

>K. 

Thus, (17) holds for ^ : = min{is:i,A:} and clearly K>0 for all 
A>0 sufficiency small. Thus, the PRFM satisfies (BR). Pick an 
arbitrary T>0. Applying Lemma 1 implies that there exists 
>7i =t]i{T)>0 such that 



Ci =£, 
Ci 

c„ = 



i = 2,...,n-l, 

li-i J 2 (20) 



1 

qn-i 



1 ]L„ 



and there exists a = a(£,ri ,^2) < 0 such that c, <a for i = 2, . . . ,n. 
2.ifD = diag(c/i , . . . ,dn) is a non-mgative diagonal matrix, with d\ > 0, 
then 

tii^{L—D)< max{ — (ii/2,a} 



x(t,ti,a)>rj^l„, foraUf>fi+T. 

Define 7,(0 : = 1 —x„ + i _,(/), ;'=!,...,«. It is straightforward to 
verify that the dynamics of the j'-system is just that of the PRFM, 
up to a reordering of the rates A{t),Ai(t). This implies that the y- 
system also satisfies property (BR), so there exists /72(t)>0 such 
that 



<0 



for allee(0,di/2). 



Proof of Theorem 4. Since L is tridiagonal with zero column sums, 
we can write 



y(t,ti,a)>ri2l„, for all f > ?i + T. 

Thus, (13) holds for ^ : = min{j7i,jj2}. This completes the proof of 
Proposition 1. 

Remark 2 Note that the proof above shows that the 
requirement that the rates are uniformly separated from zero by 
di>0 (see (4)) cannot be omitted. Indeed, if di=0 then K is no 
longer positive. In fact if we allow the rates to vanish identically 
then Proposition 1 does not hold. For example, suppose that 
A„(t) = 0. Then 1„ becomes an equilibrium point of the RFM, and 
so (13) does not hold for a = 1„. 

The next subsection includes a result on diagonal scaling of a 
tridiagonal matrix. This will be used to prove that the PRFM is a 
contraction on Qj, Ce(0,l/2). 

Diagonal Scaling of a Tridiagonal Matrix 

Let Pe]^"'*" be an invertible matrix. Define a vector norm 
I'll /> ■ M"~*M+ by |z|[ p : = |Pz|[. The induced matrix measure 
is /Zi_p(^) : =n,{PAP-'). 

Theorem 4 Suppose that Z,eK"^" is a tridiagonal matrix with zero 
sum columns, and that there exist ri,r2 >0 such that 



n < Ljj < r2, for all i # j. 

Then for each £>0 there exist 

gi=<iM>i, !=i,...,«-i, 

such that for the matrix 

P=P{e) = diag{l,quqiq2, •■■ ,9i?2 ■■ ■?«-!) 
and the matrix measure 



(18) 



(19) 



-bi 
bi 
0 

0 
0 



with 



Therefore 



91^1 
0 

0 

0 



«2 

-a2-b2 
b2 

0 
0 



02 

qi 

-02 — b2 

92^2 



0 

-Qi—bi 

0 
0 



0<ri <ai,bi<r2. 



M ■ =PLP 



0 

03 

qi 

-03— bj 



-ctjj—i bfi^i cifi 
b„-i —a„ 



(21) 



-«K-1— On-l 

qn-i 

qn—\bn—\ 



Let Ci denote the sum of the elements in column i of M, with off 
diagonal elements taken with absolute value. Since the off diagonal 
elements of M are non-negative, 

c\ = {q\-\)b\ 



: =n,(PLP-') 

the following properties hold: 



Ci=[ 1 «, + (?,■ -1)6/, i = 2,...,n-l. 
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Pick any £ > 0 and define, recursively: 

b\ 

Then clearly (20) holds. Using (21) yields q\>\-\-s\, where 
S\ : =e/r2>0. Combining this with the definition of ^2 and (21) 
implies that there exists S2>0 such that ^2 > 1 + ^2. Proceeding in 
this fashion yields an s = s{e,r\,r2)>0 such that qi> \+s for all /. 
By (20), this implies that there exists (X = oc(£,ri,r2)<0 such that 
Ci < a for i = 2,... ,n. This completes the proof of Theorem 4. 



We can now prove Proposition 2. Pick t>0, ?i >0, and aeC. 
By Proposition 1, there exists f = ((T)e(0, 1 /2) such that (13) holds. 
This implies that the off diagonal elements of L(t^x{t)) satisfy 

(5i C < Lij(t,x(t)) < ^2 (1 - 0, for all i #7, t>ti+T. 

Recall that thejacobian of the PRFM is J{t,x(t)) = L{t,x(ty)-D{t), 
with dii{t) = l(t) >Si. Let 6: = 5i/2. By Theorem 4, there exists a 
matrix P = P{e), and a scalar oc = a(s, S\^,S2{^ — C))<0, such that 
/fi^£(/(r,x(0))< max{- (5i/2,a}<0 for all t>ti-\-T. Thus, the 
PRFM is contracting on tl^ with respect to the norm This 
completes the proof of Proposition 2. 
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